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Abstract. It has long been known that the radial density profiles of globular cluster systems (GCSs) in elliptical galaxies vary 
with the total luminosities of their host galaxies. In order to elucidate the origin of this structural non-homology in GCSs, we 
numerically investigate the structural properties of GCSs in elliptical galaxies formed from a sequence of major dissipationless 
galaxy merging. We find that the radial density profiles of GCSs in elliptical galaxies become progressively flatter as the galaxies 
experience more major merger events. The density profiles of GCSs in ellipticals are well described as power-laws with slopes 
(a gc ) ranging from -2.0 to —1.0. They are flatter than, and linearly proportional to, the slopes (a s ) of the stellar density profiles 
of their host galaxies. We also find that the GCS core radii (r c ) of the density profiles are larger in ellipticals that experienced 
more mergers. By applying a reasonable scaling relation between luminosities and sizes of galaxies to the simulation results, 
we show that a w as -0.36My - 9.2, r c as -1.85Mv, and a gc « 0.93or s , where My is the total K-band absolute magnitude of 
a galaxy. We compare these predictions with observations and discuss their physical meaning. We suggest that the origin of 
structural non-homology of GCSs in ellipticals can be understood in terms of the growth of ellipticals via major dissipationless 
galaxy merging. 
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1. Introduction 

The structural properties of globular cluster systems (GCSs) 
in galaxies have long been suggested to contain fossil infor- 
mation on the early dynamical histories of galaxies and has 
accordingly been investigated both observationally and theo- 
retically (e.g., Harris 1986, 1991). In particular, radial density 
profiles of GCSs and their correlations with physical properties 
of their host galaxies have been discussed in terms of forma- 
tion and evolution of elliptical galaxies (e.g., Harris 1986; Zepf 
& Ashman 1993; Ashman & Zepf 1998; Forbes et al. 1996, 
1997; Bekki et al. 2002; 2003). In particular, observational 
studies found that the GCSs in early-type galaxies were less 
centrally concentrated than the stellar light of the host galaxy 
(e.g., Harris & Racine 1979; Strom et al. 1981; Grillmair et al. 
1994; Forbes et al. 1996), and that the outer slope of the GCS 
density profile correlated with the host galaxy luminosity (e.g., 
Ashman & Zepf 1998; Harris 1986). 

Major galaxy merging between equal-mass disk galaxies 
has long been considered to be a promising formation scenario 
for elliptical galaxy formation (e.g., Toomre 1977). Although 
this galaxy formation scenario has been discussed in many con- 
texts (e.g., the color-magnitude relation), it has not been exten- 
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sively investigated in terms of whether it can self-consistently 
explain the observed properties of GCSs around ellipticals. So 
far, only the observed higher specific frequency (Sn) and bi- 
modal color distributions of GCSs have been discussed in terms 
of hierarchical galaxy formation (Beasley et al. 2002). It is thus 
unknown what implications the observed structural properties 
of GCSs in ellipticals has for galaxy formation. 

The purpose of this paper is to propose that the radial den- 
sity profiles of GCSs in elliptical galaxies have fossil informa- 
tion on how often ellipticals have experienced major merger 
events in their dynamical histories. By using dissipationless 
N-body simulations of major merging of galaxies, we inves- 
tigate how the structural properties of GCSs in elliptical galax- 
ies change as they grow by successive (gas-free) mergers. We 
compare our theoretical predictions with observations to dis- 
cuss the origin of the observational trends. We also briefly dis- 
cuss the role of GC accretion. In particular, we focus on the 
luminosity-dependence of the GCS slope (ar gc ) which implies a 
"non-homology" in the structural properties of GCSs. 



2. Model 

We investigate the dynamical evolution of GCSs of elliptical 
galaxies formed from sequential dissipationless major merg- 
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Fig. 1. Dependences of projected number distributions of stars 
(thin) and GCs (thick) in merger remnants (i.e., elliptical galax- 
ies) on the total number of major merger events (N m ) which an 
elliptical experienced during its formation. For clarity, the den- 
sity distributions are normalized to their central values. Thin 
dotted lines represent power-law slopes (a) of a — -2.5, -2.0, 
-1.5, and -1.0. Note that the density profiles of GCSs become 
flatter for larger N m , i.e. more mergers. 



ing based on numerical simulations carried out on the GRAPE 
board (Sugimoto et al. 1990) in GRAPE 5 systems. We here 
stress that the roles of minor merging/accretion of GC in shap- 
ing radial density profiles of GCSs in galaxies are not exten- 
sively investigated in the present models. We adopt the fol- 
lowing merging scheme to model elliptical galaxy formation 
through sequential major equal-mass disk mergers. First an el- 
liptical galaxy is formed by major merging between two equal- 
mass disk galaxies. We refer to the remnant elliptical galaxy 
as the lst-generation elliptical. This 1 st-generation elliptical 
then merges with other lst-generation one to form the 2nd- 
generation elliptical. This sequential major merging can repeat 
for N m . dx times, and the structural properties of stellar compo- 
nents and GCSs are investigated for each ;'th generation ellip- 
tical (;' < N m . dx ). The total particle number used is 27670 for 
the simulations of the lst-generation elliptical and 221360 for 
those of the 4th-generation one (i.e., N mdx of 4). The initial to- 
tal particle number for GCs is set to be 1000, because an order 
of 1000 (not 10 or 100) is necessary to make robust predictions 



for the radial profiles of structure and kinematics of GCSs in 
galaxies (Bekki et al. 2005). 

Since the numerical methods and techniques we employ 
for modeling dynamical evolution of galaxy mergers with GCs 
have already been described in detail elsewhere (Bekki et al. 
2002, 2005), we give only a brief review here. The progenitor 
disk galaxies that take part in a merger (for the lst-generation 
elliptical) are assumed to have a dark halo, a bulge, a stellar 
halo, a thin exponential disk, and a GCS. The total disk mass 
and size are Md and Rd, respectively. Henceforth, all masses 
are measured in units of Md and distances in units of Rd, un- 
less otherwise specified. We adopt the density distribution of 
the NFW halo (Navarro, Frenk & White 1996) suggested from 
CDM simulations: 



P(r) 



PQ 



(r/r s )(l+r/r s ) 2 ' 



(1) 



where r, pa, and r s are the spherical radius, the central density 
of a dark halo, and the scale length of the halo, respectively. 
The dark matter to disk mass ratio is fixed at 9 for all models. 
The value of r s (typically ~ 3Rd) is chosen such that the rotation 
curve of the disk is reasonably consistent with observations for 
a given bulge mass. We adopt the R 1 ^ profile for the bulge 
with the mass of 0.17 and a scale length of Q.QARd- The radial 
(R) and vertical (Z) density profiles of the disk are assumed to 
be proportional to exp(-R/Ro) with scale length Rq = 0.2Rd 
and to sech 2 (Z/Z ) with scale length Z = 0.04/? d in our units, 
respectively. 

The GCSs in the progenitor spirals have a spherical distri- 
bution with a density profile of p(r) oc r a * c > . We adopt ar gCj i = 
-3.5 for most models (however the dependence of the final 
radial profiles of GCSs in mergers on a gCj i is also investi- 
gated). The adopted density profile of p(r) oc r~ 3 5 is consis- 
tent with that observed for the Galactic GCS (Djorgovski & 
Meylan 1994). The GCS of M31 reveals a similar surface den- 
sity distribution to the Milky Way GCS (e.g. Battistini et al. 
1993). Very little is known about the density profile of GCSs 
in other spirals. The velocity dispersion of a GCS is assumed 
to be isotropic which is consistent with the observations of the 
Galactic GCS (e.g., Freeman 1993). 

The simulations are dissipationless so no new GCs form in 
the merging process. Thus our simulations are more appropri- 
ate for relatively gas-free mergers which may be expected to 
occur at some late epochs. Recently, observational studies by 
Bell et al. (2005) have demonstrated that present day spheroidal 
galaxies with My < -20.5 on average have undergone between 
0.5 and 1 major "dry merger" (i.e., gas-poor mergers) since 
z ~ 0.7. This observational paper thus suggests that the dis- 
sipationless merger models adopted in the present paper can 
be reasonable for elliptical galaxy formation. Proto-GCs have 
been observed to form in late epoch mergers (e.g. Whitmore 
& Schweizer 1995), which implies that some fraction of GCs 
in ellipticals could be quite young. However, latest observa- 
tional studies (e.g. Cohen et al. 1998, 2003; Forbes et al. 2001; 
Kuntschner et al. 2002; Beasley et al. 2004; Strader et al. 2004; 
Pierce et al. 2005), in which more reliable age estimation of 
GCs can be done based on the higher resolution spectra and 
the improved stellar population synthesis model, have demon- 
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strated that the vast bulk of GCs in early-type galaxies are very 
old, i.e. > 10 Gyrs. These suggest that pre-existing GCs domi- 
nate in old merger remnants. 

In all of the simulations, the orbit of the two galaxies is 
set to be initially in the xy plane and the distance between the 
center of mass of the two is assumed to be 6 in our units. The 
pericenter distance (r p ) and the eccentricity (e p ) in a merger are 
assumed to be free parameters that control orbital energy and 
angular momentum of the merger. The spin of each galaxy in 
a merger is specified by two angles 0,- and <p h where suffix ;' is 
used to identify each galaxy. Here 0,- is the angle between the 
z axis and the vector of the angular momentum of a galaxy, 
and <pi is the azimuthal angle measured from the x axis to the 
projection of the angular momentum vector of a galaxy onto 
the xy plane. We specifically show four different and represen- 
tative models of merger sequences with N mdx = 4 in which 
merger pairs take one of the following orbital configurations 
(i.e., galaxy inclinations with respect to the orbital plane): A 
prograde-prograde model represented by "PP" with 6\ = 0, 02 
= 30, (pi = 0, (p 2 = 0, r p = 1.0, and e p = 0.72 a retrograde- 
retrograde ("RR") with Gi = 180, 2 = 150, 0i = 0, (p 2 = 0, r p = 
1.0, and e p = 0.72 and a highly inclined model ("HI") with 6\ 
= 30, 2 = 120, (p x = 90, and 2 = 180, r p = 1.0, e p = 1.0 and 
the low orbital angular momentum model ("LA") with 0i = 30, 
02 = 120, (p x = 90, and 2 = 180, r p = 0.2, e p = 1.0. 

In analyzing the projected radial density profiles of GCs 
(E gc ) and stars (S s ), we assume that they are approximated as 
E gc oc R a & and Z s oc R a \ where R is the projected distance from 
the center of a galaxy. We derive the slopes of a gc and a s and 
the difference of the two a gc - a s , for 5 < 7?(kpc) < 20 so that 
the derived values can be directly compared with observations. 
We determine a core radius (r c ) for the GCS by adopting an 
isothermal profile of the form E gc oc (r 2 + r c 2 ) which can be 
directly compared to the observations of Forbes et al. (1996). 
For an isothermal profile r c « 0.22R^ gc , where R^ gc is a half- 
number radius of the GCS. In order to estimate the total V-band 
magnitude of a merger remnant, we assume an initial disk mass 
of M A = 4.0 x 10 10 M o , an initial disk size of R d = 14 kpc and 
Md/Ly = 5.0, where Ly is the total V-band luminosity of the 
initial disk. The remnant's total mass is 2 Nm x Mj, where N m is 
total number of major merger events that the remnant elliptical 
experienced. 

All the calculations related to the above dissipationless evo- 
lution have been carried out on a GRAPE board (Sugimoto et 
al. 1990) in the GRAPE 5 system at the National Astronomical 
Observatory of Japan. The parameter of gravitational soften- 
ing for stellar particles is set to be fixed at 0.039 in our units 
(0.68 kpc). The time integration of the equation of motion is 
performed by using the 2nd-order leap-flog method. 



3. Results 

3.1. Structural non-homology 

Figure 1 plots the radial density profiles of stars and GCs in 
ellipticals variation with N m in the HI model sequence. The 
reference GC profile with a gc = -2.5 in this figure corresponds 
to the Galactic GC profile (Zinn 1985) which shows a gc ~ -2 
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Fig. 2. Comparison between the observed a gc - My relation 
and the simulated one. The observational data, represented by 
open circles, come from the table in the Appendix of Ashman 
& Zepf (1998). Solid, dotted, short-dashed, and long-dashed 
lines show the results of the HI, PP, RR, and LA models, re- 
spectively. Note that both simulations and observations show a 
tendency for flatter density profiles in brighter ellipticals. 
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Fig. 3. The dependence of the GCS profile slope (a gc ) and that 
for the galaxy stars (a s ) in 16 simulated merger remnants, i.e. 
elliptical galaxies. The dotted line indicates a. 



gc 



out to R = 10 kpc and ar gc 3 beyond R — 20 kpc (Harris 

1976). This figure shows that up to N m = 4, the ellipticals have 
positive ag C - a s (e.g. 0.58 for N m = 1), and thus show flatter 
profiles for GCs than for field stars, which is consistent qualita- 
tively with observations (e.g., Harris 1986, 1991). The reason 
for this is that GCs in a merger are more 'puffed up' relative 
to the stars after they inevitably absorb orbital angular momen- 
tum and kinetic energy during the merger, because most GCs 
are initially located in the galactic outer parts where such ab- 
sorption is most likely to occur. The radial density profile of a 
GCS in an elliptical is always flatter in the inner region (R < 20 
kpc) than in the outer one (R > 20 kpc). Again, qualitatively 
consistent with observations (e.g., Forbes et al. 1996). 

Both stellar and GC profiles in an elliptical become pro- 
gressively flatter as the elliptical experiences more major 
merger events (i.e., larger N m ). For example, in the HI model 



a gc is estimated as -1.99 for N„ 



1 and -1.23 for Mr 



4, 



which implies that the slope a gc of a GCS in an elliptical has 
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Fig. 4. Comparison between the observed r c - My and the 
simulated ones. Here r c is the core radius of a GCS's distribu- 
tion. The observational data represented by open circles come 
from Forbes et al. (1996). Solid, dotted, short-dashed, and long- 
dashed line show the results of the HI, PP, RR, and LA models, 
respectively. 



fossil information on how many times the elliptical has expe- 
rienced dissipationless major merger events in its formation 
history. The derived slopes for ellipticals are generally flatter 
than that observed for the Galactic GCS which has a gc « -2.5 
(Zinn 1985). This suggests that major merging can transform 
the initially steeper density profiles of GCSs in spirals into the 
flatter ones observed in ellipticals. The general trends shown in 
Figure 1 do not depend on model parameters. 

Figure 2 shows how a gc depends on the total V-band mag- 
nitude of ellipticals for four different sets of models (i.e., PP, 
RR, HI, and LA). Irrespective of the models, GCSs in brighter 
ellipticals show flatter density profiles (i.e., larger a gc ), which 
is reasonably consistent with observations. This is due to the 
fact that a gc becomes larger each time a major merger occurs. 
A least square fit to the simulation data shown in Figure 2 gives 
a gc « -0.36M V - 9.2, which is similar to the observed relation 
of a gc oc -0.3M V (Harris 1986). This result suggests that the 
origin of the luminosity-dependent a gc (Harris 1991; Ashman 
& Zepf 1998) can be understood in terms of the growth of el- 
lipticals via major dissipationless mergers. It should be how- 
ever noted, that the simulated range of a gc , for a given lu- 
minosity, is quite narrow so that the present models can not 
simply explain the scatter observed for a gc (in particular, for 
-22 < M v < -21). 

Figure 3 shows a positive correlation between a gc and a s 
(a gc « 0.93a s + 0.37), which means that ellipticals with flatter 
distributions of stars are highly likely to show flatter distribu- 
tions of GCSs. Figure 3 also shows that a s is always smaller 
than a gc (i.e., a gc - a s ~ 0.4) and thus confirms that ellipti- 
cal galaxies formed by major merging can have GCSs with the 
density profiles flatter than those of stars. The density profiles 



of the GCSs can be fit by a power-law with slope a gc « -0.4 
for the inner region of 0.5 < R < 10 kpc. However the shal- 
lower profiles of GCSs in the inner regions can not be regarded 
as "flat cores" (a gc ~ 0) which appear to be common in the 
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Fig. 5. Dependences of projected number distributions of GCs 
in merger remnants (i.e., elliptical galaxies) on the initial slope 
of radial density profiles of GCSs in spirals (i.e., a gc ,i) for the 
models with N m = 1 . For clarity, the density distributions are 
normalized to their central values. Thin dotted lines represent 
power-law slopes (a) of a = -2.5, -2.0, -1.5, and -1.0. Note 
that there are no remarkable differences between the three mod- 
els. 

GCSs of giant ellipticals (e.g., Harris 1986, 1991; Forbes et al. 
1996). 

Figure 4 shows a clear trend of larger GCS core radii (r c ) 
for brighter ellipticals, which is consistent qualitatively with 
the observational results of Forbes et al. (1996) for the overall 
GCS. The derived trend is due to the fact that r c becomes larger 
by a factor of ~ 1.4 after each major merger. We can estimate 
the luminosity-dependence of r c as r c oc -1.41M V for the ob- 
servational data with My < -20.5 of Forbes et al. (1996), and 
as r c oc -1.85My for the simulation data set shown in Figure 
4. The derived luminosity-dependence of r c is therefore steeper 
than the observed one. 

3.2. Robustness of the results 
3.2.1 . Dependence on a gc + 

Although we adopt the slope of the Galactic GCS as the initial 
radial profile for GCSs in our progenitor spirals (i.e., a gc j = 
-3.5), our results do not depend strongly on the initial profile 
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Fig. 6. The projected number distributions of stars (thin) and 
GCs (thick) in merger remnant (i.e., elliptical galaxies) with 
N m = 5. For clarity, the density distributions are normalized 
to their central values. Thin dotted lines represent power-law 
slopes (a) of a = -2.5, -2.0, -1.5, and -1.0. 

slope. Figure 5 shows how the radial profiles of GCSs depend 
on a gc j for the models with N m = 1 . It is clear from this figure 
that there is no remarkable difference in the profiles between 
the three models for R < 20 kpc: The difference of the final 
a gc between the models with a gc j = -3.0 and -4.0 becomes 
at most ~ 0.5, if we fit the profiles to the power-law ones for 
R < 20 kpc. It should be also stressed that the model with 
a gc> i = -3.5 and -4.0 shows a slightly steeper profile for the 
outer parts of the merger remnants (R > 40 kpc) compared 
with the model with a gc j = -3.0. This tendency can be seen in 
other models with different orbital configurations. 

The derived weak dependence on a gc j is due to the fact that 
violent relaxation during major galaxy merging can effectively 
wipe out the original density distributions of GCSs for -4.0 < 
Cgci ^ -3.0. The resulting profiles of a gc ~ -2 of GCSs in 
merger remnants with N m = 1 follow their dark matter halos, 
and furthermore do not depend on initial orbital configurations 
of galaxy merging. These results, combined with those shown 
in Figure 1, imply that the present results on the dependence of 
a gc on N m (thus on My) does not depend strongly on a gc ; and 
thus can be regarded as fairly robust. 

3.2.2. N m >4 

It is of interest to investigate whether the derived dependence 
of cK gc j on N m can be seen in models with N m - dx > 4 and thereby 
confirm that results are true for a possible range of N m . We note 
that the most massive ellipticals, with My < -23 mag, are typ- 
ically cD galaxies at the centres of clusters, for which accretion 
of GCs from cluster member galaxies are important (e.g., Bekki 
et al. 2003) and may well play a role in determining the radial 
profiles of GCSs. The luminosity range (for N m - dx < 4) in the 
present simulations is thus appropriate for most ellipticals but 
not cD galaxies. 

Figure 6 shows the radial density profiles of stars and GCs 
in the model with N m = 5 for the PR orbital configuration. It 
is clear from this figure and Fig. 1 that the profile of GCS be- 
comes even flatter than that in the model with N m < 4. This 
result can be seen in other models with different orbital config- 
urations. For example, a gc is estimated as -0.77 for the PR and 



-1.03 for the PP models. These results confirm that radial den- 
sity profiles of GCSs of ellipticals formed by sequential major 
merging depend strongly on N m . We note that for N mlix > 5 
and Ma = 4.0 x 1O 1O M , the resulting merger remnant has 
My < -24 mag and thus few if any observational counterparts. 

3.2.3. Spatial coverage of GCSs 

It should be noted that the spatial coverage and number of GCs 
observed for the estimation of ar gc (e.g., Ashman & Zepf 1998) 
varies from galaxy to galaxy. The spatial coverage of the sim- 
ulation data, on the other hand, is fixed at 5 < R < 20kpc for 
the a gc estimation, which suggests that the present compassion 
is not fully self-consistent. It is accordingly important for the 
present study to show how the a gc values depend on the spatial 
coverage. It is found that a gc is « -2.0 for 5 < R < 20kpc, 
« -1.6 for 5 < R < lOkpc, and * -2.4 for 10 < R < 20kpc 
in the HI model with N m = 1. This result suggests that (1) 
or gc can be larger (i.e., GCS density profiles are flatter) if a gc 
is derived for the inner regions of ellipticals and (2) a gc can 
be smaller (i.e., GCS density profiles are steeper) if a gc is de- 
rived for the outer parts of ellipticals. Given the fact that these 
tendencies can be seen in other models (PP, RR, and LA) the 
present results imply that a more careful comparison between 
observations and simulations is necessary in deriving physical 
meaning of the GCS density profiles in ellipticals. 

4. Discussion 

In the present study we suggest that sequential major merging 
has the physical effect of flattening the radial density profiles 
of GCSs so that the radial profiles of GCSs in ellipticals con- 
tain fossil information about their past merger history. Neither 
of the above two points has been previously suggested by the- 
oretical studies. 

The GCSs in ellipticals are thought to be subject to dy- 
namical destruction processes (e.g. Baumgardt 1998; Fall & 
Zhang 2001; Vesperini et al. 2003), which are not included in 
this work. However there is currently little convincing observa- 
tional evidence for destruction of GCs in ellipticals, e.g. Harris 
et al. (1998) found no variation of the bright end of the GCLF 
with galactocentric radius in M87. Vesperini et al. (2003) also 
found no variation in the GCLF with radius in M87, which they 
claimed was consistent with dynamical destruction if the GC 
mass function initially had a bell-shaped distribution. 

To date, theoretical studies which focus on GC dynami- 
cal destruction have not yet investigated the luminosity depen- 
dence of the GCS density slope (a gc ). Thus it remains to be 
seen whether GC destruction can account for this dependence, 
which is modelled here by sequential merging. 

Destruction processes are expected to be stronger in the 
galaxy inner regions and thus play some role in determining 
the inner core radius of the GCS. However it is difficult to un- 
derstand how such processes could explain the trend of r e with 
galaxy luminosity given that destruction processes are expected 
to be less efficient in the central regions of more massive ellip- 
ticals (Murali & Weinberg 1997). 
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Although the present sequential merger model can qual- 
itatively explain some of the observed GCS structural prop- 
erties, it has the following three disadvantages: Firstly, it re- 
mains unclear why GCSs in some intermediate luminosity el- 
lipticals (M v < -21.5) can show steeper (a gc < -2.0) profiles, 
which are not fully consistent with the observation (See Fig. 2). 
Secondly, the simulated core radius (r c ) of a GCS for a given 
luminosity is appreciably larger than the observed one (Note 
that the derived luminosity-dependence of r c is steeper than the 
observed one in Figure 4). Thirdly, although the derived trend 
that a gc is always larger than a s is broadly consistent with the 
observed trend (e.g., Harris 1991), GCSs in some ellipticals 
show ff gc w ff s (e.g., Ashman & Zepf 1999), which is not fully 
consistent with the present results. 

What additional physics is required to overcome these lim- 
itations? As mentioned above, the inclusion of destruction pro- 
cesses may be important but this is unlikely to explain the 
trends with host galaxy luminosity. The present study is dis- 
sipationless and so does not include the formation of the new 
GCs, which would be formed preferentially in the central re- 
gions of gaseous galaxy mergers along with field star formation 
(e.g. Bekki et al. 2002). A future study should investigate the 
role of new GC formation in shaping the radial density profiles 
of GCSs. Given our lack of understanding of the physics of 
GC formation, this is a difficult task. As we have noted above, 
the contribution of new GCs to the existing GCSs of ellipticals 
appears to be small. 

Bekki et al. (2003) showed that the strong tidal field in 
galaxy clusters can strip GCs from cluster ellipticals and con- 
sequently steepen the radial density profiles of the remaining 
GCS in the donor galaxy. Consequently, the recipient galaxy 
may obtain a flatter GCS profile. 

Recent theoretical and numerical studies suggest that if cos- 
mic reionization can strongly suppress the formation of old, 
metal-poor GCs in dwarfs at high redshifts (z > 6), then the 
radial profiles of GCSs in galaxies would be influenced at very 
early epochs (Santos 2003; Bekki 2005). 

Finally, the accretion of GCs via minor mergers (discussed 
in the next section) may also have an important effect on some 
galaxies. 

Thus the radial density profiles of GCSs in galaxies could 
be influenced by several physical processes associated with 
galaxy formation and evolution. Future studies should there- 
fore seek to determine the relative importance of each physical 
mechanism in shaping the radial density profiles of GCSs. For 
this purpose, it is vital that future observational studies pro- 
vide observational data on structural, kinematical, and chem- 
ical properties of GCs which can be compared with any the- 
oretical predictions. For example, systematic observations of 
the radial dependence of GC luminosity functions in ellipticals 
over a range of galaxy luminosities will provide a clue as to the 
dominant physical processes. 

4.1. Minor Mergers 

Although minor galaxy merging between a spiral and a dwarf 
can not form an elliptical (e.g., Walker et al. 1996; Bekki 1998), 
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Fig. 7. Final distributions of stars (left) and GCs (right) pro- 
jected onto the x-y plane for the disk galaxy model after minor 
merging of the disk with a dwarf elliptical galaxy (dE) with 
Ms = -16 mag. GCs in the right panel for this minor merger 
model are those stripped from the dE and represented by filled 
squares. The center of the disk galaxy with the initial size of 
17.5 kpc, which is not shown within the two frames for clarity, 
is coincident with the center of each frame. Note that most GCs 
are located in the outer part of the disk. 
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Fig. 8. The projected number distributions and GCs in the ini- 
tial disk model (left) and the minor merger one (right) where a 
dwarf elliptical galaxy (dE) merges with the disk. For clarity, 
the density distributions are normalized to their central values. 
Thin dotted lines represent power-law slopes (a) of a = -2.5, 
-2.0, -1.5, and -1.0. Note that the outer profile (R ~ 20 kpc) 
of the GCS of the disk becomes flattened after GC accretion 
associated with the minor merging. 



the accretion of GCs from the dwarf onto the spiral (e.g. Forbes, 
Strader & Brodie 2004) can possibly change the original den- 
sity profiles of the GCS of the spiral. Since full discussions 
on this matter are obviously beyond the scope of this paper, 
we here briefly discuss how GC accretion events can change 
the radial profiles of GCSs of galaxies in a more general way. 
Although the following discussions are on GC accretion onto 
disk galaxies, the basic results can be also applied to GC accre- 
tion onto elliptical galaxies (Bekki & Forbes 2005). Here we 
model a dwarf elliptical (dE) composed of dark matter, stars, 
and GCs and investigate the dynamical evolution of GCs during 
its orbital evolution around the present disk galaxy model. The 
full details of the dE model orbiting a spiral like the Galaxy or 
M31 can be found in Bekki & Freeman (2003), Bekki & Chiba 
(2004), and Bekki & Forbes (2005). 
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Figure 7 shows the final distributions of stars and GCs ini- 
tially within a dE with R e = 0.8 kpc, Mb = -16 mag, and 
5n = 5 (corresponding to 20 GCs in the dE), where S N is the 
specific frequency of GCs. In this model, the dE is assumed to 
be initially located 35 kpc from the center of the disk galaxy 
and have a circular velocity of 209 km s _1 at this position. The 
dE can be almost completely destroyed within ~ 6 Gyr of its 
evolution around the disk galaxy, and GCs initially within the 
dE are all tidally stripped and dispersed into the halo region of 
the disk galaxy. The tidal destruction of the dE can happen well 
outside the disk in this model so that most GCs of the dE can 
be distributed in the outer part of the disk's halo. Consequently, 
the radial density profile of the GCSs composed both of the 
original GCs of the disk and of the stripped GCs from the dE 
in the disk can become flattened, as shown in Figure 8. 

These results imply that GC accretion associated with de- 
struction of dwarfs in disk galaxies can flatten the radial pro- 
files of GCSs in the disk galaxies. These results furthermore 
suggest that although minor merging can not transform disk 
galaxies into ellipticals (Walker et al. 1996; Bekki 1998), the 
radial profiles of GCSs in disk galaxies can be significantly 
changed by minor merger events if such events are repeated 
in the dynamical histories of the galaxies. The density profiles 
of dark matter halos of galaxies can be the key parameters in 
the flattening of the density profiles of GCSs by GC accretion 
(Bekki & Forbes 2005, in preparation). Therefore the flattening 
processes by GC accretion are unlikely to depend on morpho- 
logical types (e.g., Sp, E, and cD) of luminous components, 
if all galaxies have similar dark matter profiles like the NFW 
profile. Recent numerical simulations have suggested that GC 
accretion from cluster member galaxies onto the central giant 
cDs in clusters of galaxies is highly likely (e.g., Bekki et al. 
2003). The results shown in figures 7 and 8 therefore can sug- 
gest that GC accretion can be also responsible for the origin of 
the very flat density profiles of GCSs in cDs. 



5. Conclusions 

We have demonstrated that the observed luminosity depen- 
dence of the GCS slope (a gc ) can be reasonably well repro- 
duced in elliptical galaxies formed through sequential dissipa- 
tionless major merger events. Previous numerical simulations 
have already demonstrated that the origin of structural non- 
homology in elliptical galaxies can be closely associated with 
the dynamics of major galaxy merging (e.g., Capelato et al. 
1995; Bekki 1998; Dantas et al. 2002). We suggest that struc- 
tural non-homology seen in both the stellar component and 
GCSs of ellipticals can have a common origin, i.e. the growth 
of elliptical galaxies through dissipationless major merging. 
The derived My - r c and a gc - a s relationships are also the 
direct result of this growth. Thus the radial density profiles of 
GCSs in elliptical galaxies can be regarded as containing fos- 
sil records of their merging histories. Given the fact that dy- 
namical non-homology can be closely associated with the ori- 
gin of the Fundamental Plane (e.g., Djorgovski & Davis 1987) 
of elliptical galaxies, it is an interesting observational question 
whether GCSs show an analogous "Fundamental Plane". 
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